Large violation of Wiedemann Franz law in Luttinger liquids 
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We show that in weakly disordered Luttinger liquids close to a commensurate filling the ratio 
of thermal conductivity k and electrical conductivity a can deviate strongly from the Wiedemann 
Franz (WF) law valid for Fermi liquids scattering from impurities. In the regime where the Umklapp 
scattering rate IV is much larger than the impurity scattering rate Fj mp , the Lorenz number L = 
k/((tT) rapidly changes from very large values, L ~ rjy/ri mp 3> 1 at the commensurate point to 
very small values, L ~ I^mp/TV <C 1 for a slightly doped system. This surprising behavior is a 
consequence of approximate symmetries existing even in the presence of strong Umklapp scattering. 

PACS numbers: 71.10.Pm,72.15.Eb,72.10.Bs,73.50.Lw 



In a Fermi liquid, a quasi particle carries charge e and 
has an energy of the order of fc/jT. These basic properties 
are reflected in the Wiedemann-Franz (WF) law [l], Q : 
the ratio of the thermal conductivity divided by the tem- 
perature T and the electrical conductivity, the so-called 
Lorenz number, 
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takes a universal value Lq. The WF law, L — Lq, is 
valid and routinely observed in the low-T regime of Fermi 
liquids where impurity scattering dominates. 

Deviations from the WF law, L/Lq ^ 1, in the low- 
T regime, which have e.g. been reported for high- 
temperature superconductors Q or close to quantum- 
critical points [J] , are regarded as evidence that the low- 
energy excitations cannot be viewed as electronic quasi 
particles. But even if a description of thermal and elec- 
tric transport in terms of Fermi liquid quasiparticles is 
possible, the WF law will not be valid if inelastic scatter- 
ing processes dominate which in general relax heat- and 
charge currents differently. Typically, these corrections 
to L/Lq are of the order of 1 and not very large [1, @|. 

Large violations of the WF law usually reflect a dra- 
matic change of the excitation spectrum associated with 
the opening of a gap. For example, in a Mott insulator 
a is exponentially small while heat can still efficiently be 
transported by spin fluctuations. The opposite case oc- 
curs in a superconductor where a = oo while n remains 
finite at finite T due to thermally excited quasi particles. 

In this paper, we show that small changes in the dop- 
ing can trigger enormous changes of the Lorenz number 
L in Luttinger liquids in situations where the Umklapp 
scattering rate Yjj is larger than the impurity scatter- 
ing rate, r; mp <C IV, see Fig. [T] This happens even 
in regimes where Umklapp scattering does not open a 
charge gap. This peculiar behavior can be traced back 
to the presence of approximate symmetries of the clean 
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FIG. 1: Lorenz number, I/Io, as a function of doping 8v 
away from 1/3 filling (5 = 3v c G8v/(nT)), using the variables 
of Eq. J10) (for K c = 0.6, K s = 0.8, v s /v c = 0.5). If disorder 
dominates, D > 1, L/Lq is of order one and doping indepen- 
dent. For a clean system with D <C 1, the WF law is strongly 
violated. A pronounced peak of height 1/D and width \ D 
at the commensurate filling is followed by a pronounced mini- 
mum. Inset: 8 dependence of k,q/(TLo<jo(T)), k/(TLo&o(T)) 
and a/a {T) for D = 0.001, <j (T) = (vla 2ri "- z ) g 2 ){v c j aTf 
with (3 = K s n 2 s + K c n 2 c - 3. 



system which affect charge- and heat current in a com- 
pletely different way. This has to be contrasted with a 
situation where impurity scattering provides the dom- 
inant relaxation mechanism for both heat- and charge 
currents. For this case Li and Orignac J[ have shown 
that only violations of order 1 of the WF law exist. 

When investigating the thermal or electrical conductiv- 
ity of low-dimensional systems, it is important to account 
for the role of symmetries and conservation laws even if 
these are only approximate. For example in integrable 
one-dimensional models, conductivities are usually infi- 
nite at finite T 0] as the conservation laws protects the 
currents from decaying. Small perturbations render the 
conductivity finite, but still large [§|. Below we demon- 
strate the implications on the thermoelectric effects. 
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We consider a weakly disordered one-dimensional (ID) 
metal described by a single band with the filling v = 
vq + 5v, and the electron density 2v, where vq = m c /n c 
with integers m c , n c is a commensurate filling. The low- 
energy Hamiltonian is given 0] by 



H — HlL + Hu + -ffimp 



(2) 



(27ra) 



— / dxr)(x ) (V^^ cos (-y/2> s ( a;) ] + /i.e. ) 



where i?LL is the usual Luttinger liquid Hamiltonian 
expressed in terms of spin (s) and charge (c) den- 
sities d x (f> c , s and their conjugate variable d x 9 CtS with 
[0 CjS (x), d x i9 CjS (x')] — i-k8{x — x'). Hjj is the dominant 
Umklapp scattering process where Afc = 2n c kF~m c G — 
n c G8v (with G — — ) is proportional to the deviation 
from commensurate filling and n s — 0, 1 for even and 
odd n c , respectively. The term H lmp with a Gaussian 
correlated impurity potential, (rj(x)r](x')} = D5(x — x'), 
describes a weak backscattering due to disorder. 

Even in the presence of Umklapp scattering, an ap- 
proximate symmetry closely related to momentum con- 
servation exists [loj . The so-called pseudo momentum 

P = Pt~^(N R -N L ) = P+^(N R -N L ) (3) 
Zn c 2n c 

(where N R ( L j is the number of right(left) movers) com- 
mutes with Hll + Hu (even if effects like band cur- 
vature or a weak three-dimensional coupling are added 
[lol . 11 1). Here Pt is the crystal momentum and P = 
Pt — kp{Nfi — Nl) measures the momentum relative to 
the two Fermi points. 

Because of the pseudo momentum conservation, even a 
strong Umklapp scattering may not be sufficient to relax 
the heat and charge currents. To capture this, one needs 
a transport theory which properly accounts for the role of 
conservation laws and the associated vertex corrections. 
For the non-linear interaction describing Umklapp scat- 
tering in Luttinger liquids the memory matrix approach 
to transport [12j is to our knowledge the only available 
method, especially as there are presently no numerical 
methods to calculate conductivities at finite but low T. 
As discussed in Ref. []]| , this method allows to calculate 
lower bounds to a and n in the perturbative regime, and 
gives precise results as long as the relevant slow modes 
are included in the calculation. It was shown to capture 
prominent features of observable transport phenomena, 
e.g. magnetothermal transport in spin-chains [3]. 

The first step to set up the memory matrix for- 
malism, is to list a number of relevant operators Ji 
which in our case includes the electrical current J\ = 



J c = v c K c (Nr — Nl), the heat current J 2 = Jh — 
~ Si= c s I v ^d x 4>id x 9i and the momentum operator J3 = 
P = - J2i=c, s I d x (j>id x 9 l . To leading order in Hu, H imp , 
the matrix of conductivities is then obtained from 



with the 3x3 memory matrix M = Mjj + M; mp . As 
the time derivatives dtJi = i[H,Ji] are already linear 
in the weak perturbations gu and 77, the correlators are 
evaluated with respect to All- X is the matrix of static 
susceptibilities Xij — {Ji\ Jj)uj=o with 
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The Umklapp contribution to Eq. (|4|) is given by 

2v 2 n 2 K 2 F 00 v c n c K c F- 



M, 



U 



cu?u 



where cjj 
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(6) 
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W 2n °- 

K a n 2 +K s n 2 ~l 



and = 



F mn are the dimensionless 



functions 
F„ 



dxdt te i5x (d™f c (x,t))(d2f s (x,t)) (7) 



f c (x,t) = (sinh(a; + it) sinh(a; — it)) 



^3 



(sinh (xv c /v s + it) sinh (xv c /v s — it)) 
+ {v s /v c ) 2 F 02 + (1 + (v s /v c f)F, 
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Fi = 7r[F 20 + (v s /v c ) 4 F 02 +2{v s /v c ) 2 F 1 
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which depend on doping and T via 5 = v c Ak/(jrT). 
Note that Mjj has a vanishing eigenvalue reflecting that 
[Hu,P] = 0. The disorder contribution is given by 



M, 



imp 
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nally, cr, n and L of Eq. JTJ) are obtained from 



imp 



ts , ts j r> (K a v 2 +K s v 2 )K t x^- 

K c + K s and K = nlJ\ ■ Fi- 

s v c v B (l+Kt) 



Cll, 



kq — TS a 



T 



022 



CTll 



(9) 



It should be noted that k is measured experimentally in a 
setup where the charge current vanishes, resulting in the 
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thermoelectric counter terms of Eq. ([9]). S — 621 /(Tan) 
is the thermopower. 

For given Luttinger liquid parameters v c ^ s ,K cs , the 
Lorenz number depends only on two dimensionless quan- 
tities, describing the ratio of renormalizcd disorder 
strength and Umklapp scattering and the doping: 



r- 

jj _ L imp 



Da 271 '' 3 



Tu g 2 {aT/v c )^ 



S = 



v c Ak 
nT 



(10) 



with 7 = {n\ - 1)K C + (n 2 - l)K s - 1. Fig. [T] shows 
the striking doping dependence of c, k and the Lorenz 
number L/Lq for the filling 1/3 (n c = 3, n s = 1). For 
large effective disorder, D > 1, L/Lo is of order 1 and 
there is essentially no doping dependence. For fl < 1 
one obtains instead a huge and sharp peak of height 1/ D 

and width \[T) followed by a wider dip located at i5 ~ 1 , 
where the minimum scales as D. 

This behavior can be understood by investigating the 
relation of the currents Jh and J c to the approximately 
conserved P, Eq. ((3|). From the continuity equation, 
one can show ll| that the cross susceptibility of J c and 
P is (up to exponentially small corrections) given by the 
doping 8v away from the commensurable point 



25i> 



2n c K c v c 



X31 



(11) 



while Xj h p ~ T 2 > 0. Xjp measures the "overlap" of 
the current and the conserved operator. A vanishing \ 
implies that the operators are orthogonal to each other, 
i.e. the current is not protected by the conservation law 
and can decay rapidly by Umklapp processes. Therefore, 
at the commensurate point where x.j p = 0, J c can decay 
by Umklapp processes, while Jh is protected. Indeed, as 
shown in the inset of Fig. [TJ at 6 = one obtains a ~ 
small, but k ~ l/r imp , resulting in L/L ~ 1/-D in 
the clean limit, Del. 

For finite doping, Xj p = & v > and therefore 
a ~ (<5t / ) 2 /ri mp grows rapidly until it becomes of the 
same order as the heat conductivity in the absence of elec- 
trothermal correction, Kq/T. In this regime, the leading 
contribution to k/T, however, of order 1/Ti mp is exactly 
canceled by the thermoelectric counter terms in Eq. ([5]) . 
The physical origin of this cancelation is that k is mea- 
sured under the boundary condition J c = 0. As the com- 
ponent of J c perpendicular to P decays rapidly by Umk- 
lapp, J c and P become almost parallel for small D imply- 
ing that effectively the heat conductivity measurement 
is performed under the boundary condition of vanishing 
P. Therefore K becomes of order 1/Tjj, and L/L ~ D. 
For neutral liquids a related effect is well known: while 
mass currents do not decay due to momentum conserva- 
tion, the heat conductivity measured under the boundary 
condition of vanishing mass currents remains finite (this 
situation is more transparent as momentum and mass 
current are proportional to each other while this is not 
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FIG. 2: T dependence of the Lorenz number for various 
dopings close to 1/3 filling using (|12|l (parameters as in 
Fig. H}. At low T disorder always dominates resulting in a 
T- independent L/Lq of order 1. At the commensurate point 
L/L ~ 1/73. Inset: k {T)/(TL v (D)), k(T) / {T L a (D)) 
and a(T)/a (D) for S = (red) and 5 = 10 (blue). Here 
a {D) = [Da 2n "- 3 /g 2 ] a v^/D with a = (2 - K c - K s )/y. 



the case for J c and P). Finally, for 5 3> 1 the Umklapp 
scattering is exponentially suppressed, both a and k/T 
are of order l/r; mp , and L/L ~ 1 [5[. 

In Fig. [5] the T dependence of the WF ratio, a and k 
are shown using the appropriate dimensionless variables 



S T _ v c 



Da 



2n c -3\ l H 



(12) 



Upon lowering T, the disorder close to 1 /3 filling becomes 
more and more important, D grows and L/L becomes 
of order 1 for low T . As explained above, for vanishing 
doping 5 = 0, a is much smaller than k/T as long as 
Umklapp scattering dominates. For finite doping, Umk- 
lapp scattering is exponentially suppressed at low T (see 
inset of Fig. [2]). However when it sets in (T > 1), it leads 
to a larger suppression of k/T compared to a due to the 
partial cancellations from thermoelectric corrections. 

While the theoretical analysis of the problem described 
above is most transparent for the filling close to 1/3, it 
is useful to study a case with direct experimental re- 
alizations. One possible candidate is the quarter-filled 
quasi-lD Bechgaard salt (TMTSF) 2 PF 6 [l5| where the 
anisotropy of the kinetic energy (t a :%:t c = 250 : 20 : 1 
meV) allows a Luttinger liquid description for T > \QQK. 
Two extra complications arise at quarter filling: first, in 
the absence of disorder the effective low-energy model, 
Hll + Hjj becomes the integrable sine-Gordon model, 
which formally has an infinite number of conservation 
laws on top of the pseudo momentum. For an analy- 
sis of transport one has to identify the leading correc- 
tions which break integrablity (see Ref. Q). Second, 
for Hll + Hjj there is a strict separation of charge 
and spin degrees of freedom the latter being not af- 
fected by Umklapp scattering. We therefore have to take 
band-curvature [16] into account, which couples spin and 
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FIG. 3: Lorenz number (lower curves), re and a (upper curves) 
for a system close to 1/4-filling, (ctpit = 10T -0 ' 56 is the fit to 
a) where K c = 0.22 (chosen to be compatible with Ref. [HI), 
K 3 = 0.8,v 3 /v c = 1/2, g = 0.1« c a™ c " 3/2 , n c = 4,n s = 0, m = 
l/v c a; D is in units of g 2 /(a 2n "~ 3 ), and T (in units of « c /a) 
is in the experimentally accessible regime. 



charge and breaks integrability: 
1 



H B c = ~ - /- / + 6d x cj) s d x 8 s d x 8 c 



(13) 



Here we have added an extra T-dependent chemical po- 
tential 6» = ^ (tf c + J^ 1 ) + ^ (K s + if,- 1 )) to 
account for the T-independent particle density 1v in a 
3D crystal. To leading order in 1/m, corrections to 

X arise only for X12 = Xi\ ~ ^ (V^c + V v «) and 
Xi.3 = X3i « ^ + V« s 3 )- As both A^ fl -iV L and P 

commute with Hbc, only M 2 2 gets an extra contribution, 
Af 2 C = 12 sXtvi K c {Ks 2 + K 2 S - 2) /<Im[(4cosh 2 (x + 
it) + 2) sinh(a; + it)~ sinh(xw c /v s + it)~ 2 smh(xv c /v s — 
it)~ 2 ]. As J c — > J c + P/m, ct is given by a = j u + 
2(Ti3 /m + CT33 / 77i 2 (the corresponding correction to Jh is 
subleading and therefore omitted). 

An example for the expected doping and T dependen- 
cies is shown in Fig. [3] for a filling close to 1/4 using 
parameters consistent with existing resistivity data for 
(TMTSF) 2 PF 6 [H. Both p(T) and <r(u) in this sys- 
tem can be explained [1^1 by Umklapp scattering in a 
1/4 filled Luttinger liquid with K c w 0.22 leading to 

p ^ g 2 T 16K c -3 a _ gee p ig _ [gj along the 

chain. Other parameters like K Sl v s , m, and, most im- 
portantly, disorder strength D, are not known experimen- 
tally. The absence of any visible disorder contribution to 
p(T) in the Luttinger liquid regime, T > 100K, allows us 
to estimate crudely D <C 0.0005 in units of g 2 /a 2nc ~ 3 . 
Our results shown in Fig. [3] strongly suggest that a large 
violation of the WF law (after subtraction of the phonon 
contribution not discussed here) should be observable in 
Bechgaard salts and similar materials. 



Qualitatively, the doping dependence of L/L for 1/4 
and 1/3 filling are similar. The WF ratio L/Lq shows a 
pronounced sharp peak of height 1/D followed by a dip 
for v c Ak ~ T. T-dependencies might differ in the two 
cases due to the different T dependence of D: whether 
1/D grows or shrinks upon lowering T depends on K c 
and K s . However, the most prominent T-dependence 
arises from the fact that Umklapp scattering is effectively 
switched off at lowest T for 5v > 0, resulting in L ~ Lq. 

We expect that the strong violation of the WF law in 
regimes where Umklapp scattering is large compared to 
disorder will not only occur for the strictly ID systems 
discussed here but even if a weak inter-chain tunneling 
(as in case of Bechgaard salts) is taken into account, as 
a small modulation of the ID bands does not affect the 
structure of approximate conservation laws, see 11( . Be- 



sides the disparate behavior of k/T and a an interesting 
finding of our study is the importance of thermoelectric 
corrections for the slightly doped system. In the regime 
where L/ L gets very small due to a partial cancelation 
of k and TS 2 a, the dimensionless thermoelectric figure 
of merit, ZT = TaS 2 / 'kq, which measures the efficiency 
of a thermoelectric element for power generation or re- 
frigeration, becomes 1, a remarkably large value [l7j ]. 
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